Extended Reaction Rate Integral as Solutions of Some General Differential Equations 
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Abstract. Here an extended form of the reaction rate probability integral, in the case of nonres- 
tIh . onant thermonuclear reactions with the depleted tail and the right tail cut off, is considered. The 

reaction rate integral then can be looked upon as the inverse of the convolution of the Mellin trans- 
forms of Tsallis type statistics of nonextensive statistical mechanics and stretched exponential as 
well as that of superstatistics and stretched exponentials. The differential equations satisfied by the 
extended probability integrals are derived. The idea used is a novel one of evaluating the extended 
integrals in terms of some special functions and then by invoking the differential equations satisfied 
by these special functions. Some special cases of limiting situations are also discussed. 
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1 Introduction 

Nuclear reactions govern major aspects of the chemical evolution of the universe. A proper 
understanding of the nuclear reactions that are going on in hot cosmic plasma, and those in the 
laboratories as well, requires a sound theory of nuclear-reaction dynamics. The reaction probability 
integral is the probability per unit time that two particles, confined to a unit volume, will react with 
each other (Haubold and Mathai, 1998). In the present article we will show that one can obtain 
the reaction rate probability integral as a solution of a certain differential equation, the technique 
used is a novel one. First we evaluate the integral and represent it as a special function (Mathai, 
1993; Luke, 1969). Then we can invoke the differential equation satisfied by this special function, 
thereby establishing the differential equation for the reaction rate integral. A study of this integral 
as a solution of a differential equation is made because the behavior of physical systems is usually 
studied with the help of differential equations and hence the differential equations derived in the 
present paper may become useful one day. 
Here we will consider the integrals of the following form: 

Nonresonant case with depleted tail: (for details of the integrals see Mathai and Haubold, 1988) 

/>oo 

l[ 5) = / x^e'^-^'dx, a > 0, b > 0, 5 > 0, p > 0. (1) 
Jo 

Nonresonant case with depleted tail and high energy cut-off. 

rd 

4 A) = / x "-i e -a X -bx " dX} a > , & > 0, 6 > 0, p > 0, d < 00. (2) 
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Do these integrals, as functions of 6, satisfy some differential equations? No such differential equa- 
tions can be easily seen from the integrals. We will show that these integrals will satisfy certain 
differential equations. The main results will be stated as follows. We will establish the following 
theorems for l[ S ^ and extended l[ S p and J^g . 

Theorem 1 

A constant multiple of the reaction rate integral in (1), under the condition | = m, m = 1, 2, ■ • ■ , 
can be obtained as a solution of the differential equation, 



[(-D m+1 - -viv-^-iv- ^) (V - £)] /( : ) = 



where, z = V = z£ f(z) = ^^l[ m "\z) and 



(2^) — 
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, a > 0, & > 0, a > 0, p > 0. 



(3) 
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When S = 1, a constant multiple of the reaction rate integral in (1), under the condition i 



m, m = 1, 2, • • ■ , can be obtained as a solution of the differential equation, 

[(_!)"»+! 

where, * = 7? = /(*) - ^/i^) and 



v(v -—)••■ (v - — — -) (v - <*)] f{z) = 0, 
v m m 



(5) 
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0,-1-,.. 





, a > 0, 6 > 0, a > 0, p > 0. 



(6) 



The proof of the theorem will be given later. 



Now we will consider a wider class of integrals as extended forms of (1) and (2). Consider the 
integrals 



AS) 
'i/3 



x"- 1 ^ + a{/3 - l)x d }- — e~ bx dx, a > 0, b > 0, 8 > 0, p > 0, j3 > 1, (7) 



r-d 

7 2? = / a?" -1 ! 1 - o(l - Z^^e" 6 * "dx, a > 0, & > 0, <5 > 0, p > 0, /3 < 1. 
Jo 



(8) 



These are obtained by replacing e ax& by [l — a(l — f3)x s ] 1 -f so that when /3 — >■ 1 we have the integrals 
in (l)and (2). Thus 7^ and I^j are the extended families of integrals in /j 5 ' and respectively. 
Here we obtain the following main results for the extended forms of the reaction probability integrals. 

Theorem 2 



A constant multiple of the extended reaction rate integral in (7), under the condition | = to, m = 
1, 2, • • • , can be obtained as a solution of the differential equation, 

[(-l) m+1 z( V + -i- - -?-) - V ( V - -) ■ ■ ■ (r, - — ) („ —)] f(z) = 0, (9) 
L p - 1 rap' to m rap- 
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where, z = a{ ' 3 J} bm , V = z-^, f(z) -- 



(2tt)T- 

~>- 2 . 



ATM 



l U l,m+1 



1/3 W 

a{p - 1)6" 



(10) 



a > 0, 6 > 0, a > 0, p > 0, /3 > 1. 

Corollary 2 

When 5 = 1, a constant multiple of the extended reaction rate integral in (7), under the condition 
^ = m, m = 1, 2, • ■ • , can be obtained as a solution of the differential equation, 



[(-ir +1 z(v + - a) - v(v - 1) • • • (V ~ ^) (V a)] f(z) = 0, 



m 



(11) 



where, z 



a(P-l)b" 



p[a M T j^ ml 

17 /iM-z) and 

(2tt)T- 



(27T) 



/>[«(/* -i)]"!^)™* 



"a(/3 - l)6 m 




m m 





(12) 



a > 0, 6 > 0, a > 0, p > 0, /3 > 1. 

Theorem 3 

A constant multiple of the extended reaction rate integral in (8), under the condition | = m,m = 
1, 2, • • • , can be obtained as a solution of the differential equation, 



(13) 



where, z = , n — z4-, f(z) = pK 1 P )1 """^ I™f(z) and 

m ' ' r (^f)(2vr)-4^ 2,3 v ; 



r (H)( 2 ^) J 

T mp( \ _ V ' 



p[a(l - /?)] m "m2 
a > 0, 6 > 0, a > 0, p > 0, /3 < 1. 

Corollary 3 
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When 6 = 1 a constant multiple of the extended reaction rate integral in (8), under the condition 

(15) 



i = m, m = 1, 2, • • • , can be obtained as a solution of the differential equation 



[{-l) m z{ V --l--a)-r,( V --)..-(r ] - — ) („ - «)]/(*) = 0, 
\ 1 — n to' m 



where, z = a(1 ~g, )b " 



dz ' ' w r (^|)(27r)-4^ 2/n ; 
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a > 0, b > 0, a > 0, p > 0, /3 < 1. 

The proofs of the above theorems will be given after evaluating the integrals first. 



2 Evaluation of the integral l[ S ^ 

The following procedure is available in Mathai and Haubold (1988). For the sake of completeness a 
brief outline of the derivation will be given here. 

/>oo 

l[ s) = / z Q - 1 e- ax *- bx ~ P dx, a > 0, b > 0, S > 0, p > 0. 



Here the integrand can be taken as a product of positive integrable functions and then we can apply 
statistical distribution theory to evaluate this integral. Let x\ and x 2 be real scalar independent 
random variables having densities 

f , n _ / cix"e~ aa: i, < xi < oo, a > 0, 5 > nr . 

o, elsewhere 10 



and 



J v ; 1 0, elsewhere, v ; 



where ci and c 2 are normalizing constants. Let us transform x\ and x 2 to u = x\x 2 and u = x\. 
Then the marginal density of u is given by 

9i(u) = [ -h(v)f 2 (-)dv (19) 

Jv v v 

= Cl c 2 / v Q - 1 e - Bv4 - bv ~ P dv, where b = u", 6 > 0, p > 0. (20) 
Jo 

Let us evaluate the density through expected values or moments. 

E{u s - 1 ) = E{x s 1 - 1 )E{x s 2 - 1 ) (21) 
due to statistical independence of x\ and X2. 

p oo 

Eix^ 1 ) = Cl / x1 +s - 1 e- ax ^Ax 1 . 
Jo 

Putting y = axf and evaluating the integral as a gamma integral, one has, 

Efr'r 1 ) = T-^ r (T + t). 5R(a + ^)>0, (22) 

where 3?(.) denotes the real part of (.). 

E(x s 2 ~ 1 )=c 2 / 
Jo 



^2 1 e x ' 2 dx 2 . 



Putting y — x 2 , we get 
From (22) and (23) 



E(x^ 1 ) = Jr(- p ), n(s)>0. (23) 



,a s 



E(u s -') = — ^3-r(-)r( 7 + -), > 0, 5R(a + s) > 0. (24) 

dpas^s pod 
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Looking at the (s — \) th moment as the Mellin transform of the corresponding density and then 
taking the inverse Mellin transform we get the density of u, 



9i(u) = 

Comparing (20) and (25) 

AS) _ 
i i — 



ClC 2 1 

Spa% 2m' 



J c — zoo pOO 



(25) 



. e -ax s -bx » dx 



POO 

Jo 

1 1 />C+ZOO 

Spa s 2m J c _ i(x p 8 5 



This contour integral can be written as an ii- function (Mathai and Saxena, 1978). That is, 



(S) 



5pa~z 



'0,2 



as bp 



, a > 0, b > 0, a > 0, S > 0, p > 0. 



(26) 



(27) 



Make the transformation 4 = s\ in (26) 



A s ) 



1 1 

pa% 2m 



_ r(^)r(| + Sl )(a6f)-d Sl 



(28) 



Let us consider the special case where - = to, m = 1, 2, • • • . In physics problems p = \ and 5 
an integer. Then this assumption of - = to, m = 1, 2, • • • is meaningful at least in some physical 



problems. That is, 



1 1 rci+ioa 

7 ( mp) = + / r( mai )r(— + Sl )(ab m r si d Sl . 

pa —2m]^_ loo y mp 



pa 

This can be reduced to a G-function by using the multiplication formula for gamma functions, 
namely, 

Timz) = (2n) k ^ L m mz -^T(z)T(z + — ) ■■■T<z+ ^— -), m = l,2,--- . (29) 



to 



Then we have, 



i[ mp) = ^Ux^ r +4 °° TO^r( Sl )r( Sl +l) • ■ . r(* 1+ — )r(* 1+ — ) (ab m )-^d Sl 

pa™? 2m J Cl _ too m' m mp' 



(30) 



On evaluating (30), we get (4). (see Mathai and Haubold, 1988; Saxena, 1960); Haubold and Mathai, 
1984). 

Particular case, 8=1 

When 8 = 1, we get the reaction rate integral in (1), under the condition - = to, 
to = 1, 2, • • • as, 

/•OO 

h = / a; a_1 e- aa; - 6x " P da;, a > 0, 6 > 0, p > 0. (31) 
Jo 

On evaluating this integral, we get (6). Observe that m — 2 is a real physical situation, see for 
example, Mathai and Haubold (1988). 
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3 Evaluation of extended integrals 



We have, 



AS) 

Z l/3 



[1 + a(/3 - l)x d ]-f^c- bx "dx, a > 0, b > 0, 5 > 0, p > 0, /3 > 1. 



As /3 — >• 1, [1 — a(l — j3)x 5 ] 1 -f> becomes e axS so that we can extend the reaction rate integrals in 
(1) and (2) using the pathway parameter /3. Then here arise two cases:(i) (3 < l,(ii) (3 > 1. 



Case(i), /3 < 1 

/■(*) 

1 2(3 

Here let us take 
and 



/ x^^l - o(l - /3)a; 6 ]T^e- bx "dx, a > 0, 6 > 0, 6 > 0, /3 < 1. 
Jo 

cia:f [1 - o(l - /3)a;i 4 ]^, < an < [^by] ^ « > 0, <5 > 0, /3 < 1 



/2(s 2 ) 



elsewhere 



c 2 e x 2 1 o < x 2 < oo, p > 
0, elsewhere, 



(32) 



(33) 



where ci and C2 arc the normalizing constants. Let us consider the case where d — [ a ^_^ ] J in (8). 
Note that /i(xi) in (32) with a — and 6 = 1 is Tsallis statistics leading to nonextensive statistical 
mechanics. Also observe that the functional part of (32) for a — and 6 = 1 gives the power law as 
well. 



_d_ 1-/1(2:1) 



dx ci 



fi(xi)- 



Cl 



(34) 



Note that /2OE2) in (33) is what is known as stretched exponential in physics literature. Thus the 
extended reaction rate integral in (8) is the Mcllin convolution of Tsallis nonextensive statistics and 
stretched exponentials. The starting publication of nonextensive statistical mechanics may be seen 
from Tsallis (1988). Then proceeding as before, we get, 



AS) 
'2/3 



r 1 1 * 

= ( ^ x a - l [l-a{l-p)x 5 ]^e- bj ' 
Jo 



dx 



r ( 1 + T3«) 1 f c+i °° T (- P )nf + l)[aHl-/3)h^}- s d S 



p5[a(l -/3)] T 27r * 

r(i + iV 



r(£L 4. 1 J — L_ 



p*[a(l-0)]' 



1,2 



a'(l-/3) 7 6p 



(f+i+T^.i) 



(oi),(f,i) 



(35) 
(36) 
(37) 



a > 0, 6>0, a>0,5>0/9>0, /3<1, 5R(s) > 0, 5R(a + s) > 0. 
Make the transformation | = Si in (36) 



A*) 

'2/3 



r^ + i 1 ?) 1 r"°°r(^)r(f + Sl )[a(i -/?)^]-^d Sl 



p[a(l-/3)] 



t 2ni 



(38) 



Let us consider the case where | = m, m = 1, 2, • ■ • . Then we get (14). 
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Lemma 1: As — > 1, I^J becomes I 2 '. 
Proof: 

i r c+i °° r(f )r(f + f)(a*6*)- r(i + ^)(i - /?)-(f+f) 



lira /if? = lim 

/3->i /p 27TI 



p<5a ■ 



r(f + 1 + i + t^; 



ds. (39) 



Now apply the asymptotic formula for gamma functions, namely, for \z\ — > oo and a is bounded 
(Mathai, 1993), 

T{z + a)^ (27r)h z+a -h- z . (40) 

Apply this to the gamma ratios in (39) by taking z as — ^ and a as 1 and f + | + 1 respectively. 
Then, 



r(i + T ^)(i~/3)-(f + f) 

hm ; ; — r = 1 

^ r(f + f + i + rb) 



Hence 



lim /, 



26 



(41) 



(42) 



which establishes the result. 
Case(ii), /3 > 1 

/■OO 

J$ = / ^^[l + a(/3 - l^n^e-^dx, a > 0, 6 > 0, <5 > 0, p > 0, j3 > 1. 
./o 



/o 

Here let us take 

A (si) 

and 



cix?[l + a(/3 - l) Xl s Y~, < < oo, a > 0, <5 > 0, /3 > 1 
0, elsewhere 



c 2 c x 2, < x 2 < oo, p > 
0, elsewhere, 



(43) 



(44) 



where ci and c 2 are the normalizing constants. Observe that fi(xi) of (43) is nothing but the 
superstatistics of Beck and Cohen (2003), and the density in (44) is the stretched exponential. Hence 
(7) can be looked upon as the Mellin convolution of superstatistics and stretched exponentials. A 
large number of published articles are there on superstatistics. Then proceeding as before, we get, 



Jo 



-bx 



where K 



r(^rjMK/3-i)]T 



e ™ dx 

( a i fi~ 2 1 1 n 

(o,i),(f>i)- 

, a > 0, 6 > 0, a > 0, p > 0, S > 0, /3 > 1, SR(s) > 0, 



= K Hf 



a? (/3 - l) a 6p 



(45) 



5ft(a + s) > 0. Make the transformation | = si and let - = m, m = 1, 2, • • ■ in (45). Then we get 
(10). 



Lemma 2: As /3 — > 1, Iyj becomes l[ S \ 
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Particular case, 6=1 

When S = 1, we get the reaction rate integral under the condition - p = to, 
to = 1, 2, • • • as, 

[ i ] 

hp = I ° (1 "" > x a - 1 [l-a(l-/3)a;]T J ire- ba:_P da;, a > 0, b > 0, p > 0, £ < 1, (46) 
Jo 

on evaluating this integral we get (16), and 

f°° -i 
hp= x a - 1 [l + a(P-l)x]m c - bx "dx, a>0, 6 > 0, p > 0, /3 > 1, (47) 
Jo 

on evaluating this integral we get (12). 

Proof of Theorem 1 

For proving the theorem we will make use of the fact that the reaction probability integral 
as well as the extended reaction probability integrals, as given in (1),(2),(7),(8), under the condi- 
tion (5 = to/j, to = 1, 2, • • • can be written in terms of G-functions. Hence when this condition is 
satisfied we can invoke the properties of G-functions. It is well known that the G- function defined by 

g?;?[az::Z] - hi*®*" 6 ** ( 48 ) 

where 

, {UUm + s)} Ol"=ir(i- fflj -- a )} 

ns) " {n]= m+1 r(i - hi .)} {UU + i + -)} ' 

aj, j = 1, 2, ...,p and bj, j = 1, 2, q are complex numbers, L is a contour separating the poles of 
T(bj + s), j = 1, 2, to from those of T(l — aj — s), j = 1, 2, n, 
satisfies the following differential equation. 

[(-If ~ m ~ n z f[( V -a i + l)-f[{r,- bj)] G(z) =0, V = z± (49) 
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This equation is intuitively evident from the following facts: 

7]Z- S = Z^{Z- S ) = {-S) Z - S - 

(77 - dj + l)z~ s = (1 - dj - s)z~ s ; 
(1 - aj - s)T(l - dj -s)= T(2 - a, - s), 
(see Mathai, 1993). The depleted case of the reaction rate integral 



I? = / 
Jo 



x^e'™ - bx "dx, a > 0, b > 0, 6 > 0, p > (50) 



.— p 



can be expressed in terms of G-functions under the conditions | = m, m = 1, 2, ■ • ■ shown in (4). 
That is, 

- ' - r ab m 

(51) 



ii 1 



fn \ — — ^ T (mp) /-fm+1,0 

(2?r) 2 m2pa mp A = G 0i?n + 



TO" 



a > 0, 6 > 0, a > 0, (5 > 0, p > 0, - p =m, m = 1, 2, • • • . 
The G-function in (51), satisfies the differential equation in (3). So the left hand side of (51) also 
satisfies the differential equation (3). Similar are the proofs of the other theorems and hence deleted. 
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